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ABSTRACT 

The discovery of Ultra-Faint Dwarf (UFD) galaxies in the halo of the Milky Way extends the faint end of the 
galaxy luminosity function to a few hundred solar luminosities. This extremely low luminosity regime poses a 
significant challenge for the photometric characterization of these systems. We present a suite of simulations 
aimed at understanding how different observational choices related to the properties of a low luminosity sys- 
tem impact our ability to determine its true structural parameters such as half-light radius and central surface 
brightness. We focus on estimating half-light radii (on which mass estimates depend linearly) and find that 
these numbers can have up to 100% uncertainties when relatively shallow photometric surveys, such as SDSS, 
are used. Our simulations suggest that to recover structural parameters within 10% or better of their true values: 
(a) the ratio of the field-of-view to the half-light radius of the satellite must be greater than three, (b) the total 
number of stars, including background objects should be larger than 1000, and (c) the central to background 
stellar density ratio must be higher than 20. If one or more of these criteria are not met, the accuracy of the re- 
sulting structural parameters can be significantly compromised. In the context of future surveys such as LSST, 
the latter condition will be closely tied to our ability to remove unresolved background galaxies. Assessing the 
reliability of measured structural parameters will become increasingly critical as the next generation of deep 
wide-field surveys detects UFDs beyond the reach of current spectroscopic limits. 

Subject headings: galaxies: dwarf, fundamental parameters - Local Group, structure - photometry 
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1. INTRODUCTION 

The structural parameters of galaxies are of fundamental 
importance. Uncertainties in estimating a galaxy's physi- 
cal size, central surface brightness, or total luminosity can 
quickly propagate into other key galactic properties. For 
example, galaxy mass measurements, based on either the 
virial theorem or Jeans eq uation modeling, depend linearly on 
the ga lactic radius (e.g.. Ivan der Marel et al Jfl994t IWoif et al. 
2010). While estimating structural properties should be 
straightforward, practical considerations such as a limited 
field size or background estimation can make this a compli- 
cated problem. 

We have entered a new regime in estimating structural 
parameters with the discovery of Ultra-Faint Dwarf (UFD) 
galaxie s around the Milky W ay and M31 dWillman et alj 
" 2005albt IZucker et alj|20 06a b: Be lokurov et all 120061 120071 
2010tlMartinetal.ll2009l) . These objects are found via sta- 



tistical over-densities of resolved stars in the Sloan Dig- 
ital Sky Survey (SDSS) and other large imaging surveys 
with the faintest ones having total luminosities of merely 
M v ~ -1.5, or 300 Up dMartin et al.ll2008t iGeha et alj|2009l: 



iBelokurov et alj 120101) . The exact nature of an individual 
UFD is often debated (is it a star cluster or dark matter- 
dominated galaxy? is it in dynamical equilibrium?) and 
structural parameters are an important tool in this assessment. 
However, at the detection limits of the SDSS, their low total 
luminosity translates into a few tens to hundreds of individual 
stars observed as UFD members. In this regime, a few stars 
can significantly alter the derived size or total luminosity of a 
given object. 

Structural parameters are determined throughout the liter- 
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ature by measuring the surface brightness of a galaxy as a 
function of r adius, and then fi tting an analytic function to this 
profile (e.g. JCaon et al.l[l9 93). The surface brightness profile 
is usually well determined for galaxies based on an integrated 
light profile, or in the case of nearby luminous (My < -8) 
Local Group galaxies , by summing many thousand s of stars 
in concentric annuli ( Irwin & Hatzidimitriou 1995). In ap- 
plying these techniques to the UFD galaxies, shot-noise due 
to the low number of available stars can produce unreli- 
able results. As an alternative to surface brightness fitting, 
iMartin et alj J200 8) determined structural parameters for the 
Milky Way UFD galaxies using a Maxi mum Likelihood (ML) 
met hod, building on pre vious work by iKleyna et alJ (119981) 
and Westfa ll et alJ (12006). In this method, a fixed analytic pro- 
file is assumed and its free parameters (typically galaxy cen- 
ter, half-light radius, ellipticity, position angle and the back- 
ground level) a re fit simultaneously using all the av ailable 
stars (see also I Sand et al] 120 09: Mun oz et alj|20101) . This 
avoids the need to bin or smooth data and should provide more 
reliable estimates for systems with low number of stars. 

There are regimes, however, where even the maximum like- 
lihood method will produce unreliable structural parameters. 
In this paper, we explore the transition region between the 
extreme case where only a few member stars of a resolved 
galaxy are identified and the well measured galaxies as de- 
scribed above. A quantitative understanding of this regime 
will be of increasing importance since large surveys, such as 
SkyMapper, Pan-STARRS, DES, and LSST, are predicted to 
discover many tens to hundreds more low luminosity systems 
dTollerud et al.ll2008l) . 

Here we perform a representative set of simulations to ex- 
plore the accuracy of the ML method for structural parame- 
ter estimation in cases when the nu mber of st ars d etected in a 
galaxy is a few thousand or less. In § 12. H and § !2.2l we describe 
the grid of simulated resolved galaxies, and the ML method 
used to determine structural parameters in these simulations. 
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FIG. 1 . — Variation of the stellar background density as a function of both 
Galactic latitude and longitude at the position of the ComBer satellite over a 
1° region. These data were simulated using the Besancon Milky Way model. 
Variations in the background counts are consistent with Poisson noise (error 
bars). 

We then consider parameter optimization for the particular 
cases of the SDSS and LSST surveys(§[3]i, as well as tar- 
geted observations with the Canada-France-Hawaii-Telescope 
(CFHT, §|?J. In §|5]we explore the accuracy of the recovered 
structural parameters split into the separate effects of density 
contrast, the total number of stars, and the effects of a limited 
field-of-view (FOV). Finally, in §|6]we summarize our results 
to distill them into a set of generalized guidelines for future 
observations. 

2. METHODS 

How well can structural parameters be measured in the 
regime where a galaxy is detected as only a low number of 
resolved stars? To answer this question, we generate a series 
of simulated UFD galaxies and assess how different observa- 
tional variables impact our ability to determine their structural 
parameters. We generate a series of simulated resolved galax- 
ies, varying the total number of stars, the relative density of 
background stars and the physical size of the object relative 
to the FOV. We then determine the structural parameters of 
these systems using a maximum likelihood method and com- 
pare our recovered values back to the true inputs. 

2.1. Simulations 

The stellar density profiles of dwarf galaxies, including 
the UFDs, are typically model ed using analytic density laws 
such a s exponential, Plummer (Plummer ]191 11) or King dKingl 
119621) profiles. We construct our model galaxies using a 
Plummer profile, but note that this choice does not affect our 
conclusions provided the same density law is used to gener- 
ate both the simulations and recover the structural parameters. 
The projected Plummer density profile is: 
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where R\ji is the half-light radius, the galaxy center is at (xo, 
yo) and Eo is the central density. As we are interested in ex- 
ploring the effects of a limited FOV on the recovered struc- 
tural parameters, we measure R\n in units of FOV, L. 



TABLE 1 

Approximate value for our different parameters for three 
representative dwarf galaxies assuming sdss and lsst-like 
surveys as well as targeted observations with megacam on 
the cfht. we also list the assumed /"-band depths and fovs. 

Survey/Telescope r-band depth L R\/i ^dwarf ^T a Sq 



ComBer (My - 


= -3.8, Rin = 6.0', d 


= 44kpc) bc : 








SDSS C 


22.5 


2° 


0.05 


99 


3310 


~5 


LSST s j n g] e _ eX p 


24 


1° 


0.1 


560 


1305 


-25 


LSST C o_ ac ided 


27 


1° 


0.1 


3420 


8605 


~22 


CFHT 


24.5 


2° 


0.05 


750 


3740 


- 35 


CFHT 


24.5 


1° 


0.1 


750 


1500 


~ 35 


Segue I (My = 


-1.5, R 1/2 = 4.4',d = 


= 23kpc) c : 








SDSS C 


22.5 


2° 


0.037 


65 


4340 


~3 


LSSTjjjngie—exp 


24 


r 


0.075 


200 


945 


~ 16 


LSSTcQ-added 


27 


i° 


0.075 


1060 


3610 


~25 


CFHT 


24.5 


2° 


0.037 


260 


3245 


~ 21 


CFHT 


24.5 


1° 


0.075 


260 


1005 


~ 21 


Boo II (My = - 


-2.7, #1/2 = 4.2', d = 


42kpc) c : 








SDSS C 


22.5 


2° 


0.035 


37 


3940 


~3 


LSST s j n g] e _ eX p 


24 


1° 


0.07 


210 


960 


~20 


LSST c0 _ ac i(jed 


27 


1° 


0.07 


1280 


6020 


- 17 


CFHT 


24.5 


2° 


0.035 


290 


3275 


-25 


CFHT 


24.5 


1° 


0.07 


290 


1035 


-25 



a The number of backgro und counts has be en derived using background den- 
sity values taken from Martin et al. (2008) for SDSS and from Munoz et al. 
12010) for the CFHT/MegaCam survey. We note that in general, the latter 
values are lower than those from SDSS, despite the deeper photometry. For 
the case of an LSST-like survey we used CMD-filtered counts from the Hub- 
ble Ultra Deep Field and assumed that 90% of background galaxies can be 

removed; 

b From lMunoz et al l fMf]) . 
c From lMartin et all <2008T) . 

Any stellar survey will also include some contaminating 
sources, either foreground Milky Way stars or unresolved 
background galaxies. We place our satellites over a ho- 
mogeneous background density with Poisson variations. 
Throughout this paper, we assume Ej, = 1 over the relevant 
area, and report the central density Eo in units of the back- 
ground density. We have assessed the validity of our assump- 
tion of a homogeneous background using both the Besancon 
model at the Galactic position of Coma Berenices (ComBer) 
and Bootes II (Boo II), and CFHT imaging of the same UFD 
galaxies. We find that variations for both datasets, in bins of 
12 arcmin across a one degree patch in the simulations and 
the one degree images, are on the order of 10% and are con- 
sistent with Poisson noise (see Figure [TJ. We conclude that 
background variations over an area of one degree, relevant for 
Milky Way UFDs, are not significant and our choice of a ho- 
mogeneous background, as a first approximation, is justified. 

Given the Plummer profile above and assuming circular 
symmetry, a simulated dwarf galaxy is fully described by 
three free parameters: R\/2, Eo and Nt- The final parame- 
ter Nt is the total number of objects, including UFD stars and 
background detections. In TableQ] we list approximate values 
for a few representative dwarf galaxies as observed by differ- 
ent instrument/survey combinations. 

In Figures |2] and [3] we show a dwarf galaxy simulated un- 
der different observing conditions. Figure [2] shows one re- 
alization of a satellite resembling the UFD galaxy Bootes II 
(Boo II, My = -2.7, d = 42kpc) as it would be observed by 
the CFHT/MegaCam combination down to a magnitude limit 
of r = 24.5. The galaxy is clearly seen above the background 
level as a centrally concentrated group of resolved objects. 
In contrast, Figure [3] shows the same galaxy as observed by 
the SDSS. Due to the much shallower data, the galaxy is dif- 
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FIG. 2. — Simulated star map of an ultra-faint galaxy with parameters: 
#1/2 = 0.07, S = 25 and N T ~ 1000, similar to a Boo II-like satellite ob- 
served with the CFHT/MegaCam combination down to a magnitude limit of 
r = 24.5. 

ficult to discern above the background counts. Indeed, the 
UFD galaxies were only recently discovered because the stel- 
lar density contrast between the background and member stars 
is low- in a typical image a UFD cannot be "seen" without ad- 
ditional filtering of the background. 

We run a grid of simulations to cover most of the range of 
UFD observations presented in TableQ] We vary Rin between 
0.025 and 0.5, So from 3 to 100 and the total number of stars 
Nt from 300 to 5000. For simplicity, the dwarf galaxies are 
placed at the center of the simulated field. 

2.2. The Maximum Likelihood Method 

We estimate the structural parameters of our simulated 
galaxies using a maximum likelihood (ML) method. We as- 
sume that the stars are distributed according to a density pro- 
file, 

E(x,y, p) = £ dwarf (x,;y, p) + E b (2) 

where Sdwarf is the Plummer profile of equation Q] and T,/, is 
an arbitrary background density we choose to normalize to 
one. The density depends both on the position in the field 
(x,y) as well as the structural parameters p = (^ 1 / 2 ,Eo,xo,yo)- 
Note that, in addition to the half light radius and central sur- 
face density, we also fit for the central position of the galaxy 
(xo,yo)- We continue to assume circular symmetry for sim- 
plicity. 

The probability £, of observing star i is then given by 

c= £(*,-,y,-,p) (3) 
J n Z(x,y,p)dxdy 

where J n is the integral over the observed field and normal- 
izes the probability distribution. The probability of observing 
the full data set given a particular choice of parameters p is 
therefore 



FIG. 3. — Same as Figure [2] but for a Boo II-like satellite as observed in 
the SDSS. This realization has parameters: R[/2 = 0.035, Eo = 3 and N = 
4000. While the galaxy is not obvious by eye, it is detected as a statistical 
overdensity of star. 

The maximum likelihood estimate of the parameters is de- 
fined to be the vector p that maximizes the probability of ob- 
serving the data. The initial values of the measured parame- 
ters were chosen arbitrarily (not at the true value); changing 
the starting position made no difference t o the final results- 
Comparing this to previous w ork (M artin et alj 120081: 
ISand et al.ll2009t iMunoz et aill2010l) . we note that the proba- 
bility is invariant under multiplicative scalings, explicitly jus- 
tifying our fixing Sb = 1 . An alternative approach is to explic- 
itly fit for the background density, in which case the normal- 
ization integral is t hen the total number of stars A^r- This is the 
approach taken bv lMartin et alj (120081) . although they make a 
further simplifying assumption that the FOV is large enough 
that integrating the dwarf profile over the field is equivalent 
to integrating it to infinity. We do not make this assumption 
here, since we are also interested in cases where the FOV is 
comparable to the size of the dwarf. 

For each set of free parameters (^1/2, So and Nt) we run 
1000 realizations and generate a distribution for each of the 
recovered inputs. To assess the quality of the recovered pa- 
rameters, we report the width of these distributions. We de- 
fine the quantity a p as the half width of the central 68% of the 
distribution. While the distributions are not necessarily Gaus- 
sian, this quantity is equivalent to the one-sigma width for a 
Gaussian distribution. In addition, it is worth mentioning that 
in all cases, no appreciable bias in the recovered parameters is 
observed, that is, the observed distribution of measured values 
is roughly centered on (within at most 10% of) the true input 
number. In the analysis presented below, we focus primarily 
on the behavior of the half-light radius and therefore report 
only o-R xp values. 

3. STUDYING UFDS IN WIDE-FIELD IMAGING SURVEYS 

In studying UFD galaxies, two different scenarios are rele- 
vant: observations made using surveys, and targeted observa- 
tions with a specific telescope/instrument combination. These 
two cases present different observational constraints and we 
examine each separately. 
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FIG. 4. — Fractional uncertainty cr^n as a function of the satellite's relative 
size, R\n for two UFDs, Boo II (solid lines) and ComBer (dotted lines), "ob- 
served" with two photometric surveys, SDSS (upper lines) and LSST (lower 
line). For each line in this figure, the density contrast has been kept constant 
and therefore the total number of stars varies from high numbers on the left 
of the lines to low numbers on the right. This figure shows that when the 
density contrast is fixed, higher values of R\n will not necessarily result in 
better uncertainties because the number of stars decreases rapidly. 

All of the UFDs have so far been discovered in the SDSS, 
but this sample is far from complete due to the survey mag- 
nitude limits and sky coverage. Upcoming surveys, such 
as SkyMapper, Pan-STARRS, DES and LSS T are expected 
to fin d many tens to hundreds more satellites dTollerud et al.l 
2008). For any wide-field imaging survey two variables perti- 
nent to estimating structural parameters are fixed: the number 
of observed stars belonging to the dwarf galaxy A^warf , and the 
density contrast. Both depend on the quality and depth of the 
photometry. For a given satellite's angular size and distance, 
the central density depends directly on the number of mem- 
ber stars detected by the survey, and is set primarily by the 
depth of the photometry, and to a lesser extent, by the qual- 
ity of the data. The background density increases with the 
depth of the survey but better quality data allows for more ef- 
ficient removal of background unresolved galaxies. The only 
free choice an observer faces is the area around the satellite 
used to carry out the structural parameter determination. This 
choice will impact both the relative size of the satellite with 
respect to the FOV and the total number of stars Nt ■ Naively, 
one expects the area around the satellite to have a mild effect, 
if any, on the accuracy of the derived structure, but in the new 
regime of ultra-faint satellites it is an essential observational 
consideration. 

We focus below on two surveys, SDSS and LSST, and ex- 
amine how selecting different area sizes affects the accuracy 
of the measured half-light radii. 

3.1. SDSS 

To illustrate how our simulations can guide the data analy- 
sis for study of UFDs, we carry out mock observations for two 
example d warf galaxies: Boo II (My = -2.7 ± 0.9, / = 353°.7, 
b = 68°. 8; Mar tin et al.ll2008l) and Coma Berenices (ComBer , 
My = -3.8 ±0.6 / = 241°.9, b = 83°.6: iMufioz et al.ll2010h . 
The former is one of the smallest and faintest of the UFDs, 
while the latter is rather "typical". We use for reference the 
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FIG. 5. — Same as Figuref4] but for the fractional uncertainty in luminosity 



structural parameters of Marti n et~aT] (120081) who used SDSS 
data to carry out a ML study of all UFDs discovered before 
2008. For their structural analysis, these authors selected stars 
within 1 deg radius of each galaxy down to the SDSS magni- 
tude limit of r ~ 22.5. 

From this catalog, Boo II has a measured half-light ra- 
dius of 4.2'_ 14 and is located at a distance of 42kpc. In 
terms of the variables studied by our simulations, this results 
in S w 3, Ri/2 = 0.035 and N T = 3940. The simulations 
show that this combination of parameters for Boo II yields 
a <jR in of rou ghly 30% (similar to the uncertainty reported by 
iMartin et al J 120081) and a ol ~ 90, indicating that the mea- 
surement based on SDSS data is hig hly inaccurate. It is worth 
pointing out that lWalsh et alJ (2008), using deeper photomet- 
ric observations of Boo II by roughly one magnitude, reported 
afli/ 7 = 2.5' ±0.8 nearly 40% smaller than the 4.2' found by 
Marti n~et"aT] d2008). consistent with the uncertainty obtained 
from our simulation s. For the case of ComBer, the results of 
IMartin et alJ (l2008h give a half-light radius of 6.0' ± 0.6 and 
a distance of 44kpc, which results in So « 5, R\/2 = 0.05 and 
Nj = 3310. This combination of parameters gives a^ p ~ 20% 
and , ol ~ 50, an improved value as compared to Boo II. In 
the cas e of ComBer, the uncertainties reported by Mart in et alJ 
(2008) in their measured half-light radii are somewhat lower 
that those from our simulations. To estimate their errors, these 
authors assumed that the likelihood function is well-described 
by a Gaussian distribution near its maximum, an assumption 
that probably breaks for the fainter objects. 

As mentioned above, potential improvements in the accu- 
racy of the measured sizes can only be achieved by chang- 
ing the area around the satellite used to select UFD mem- 
ber candidates. In particular, we explore the effect of choos- 
ing a smaller region. Given that the density contrast will re- 
main constant, using a smaller FOV will bring down the total 
number of stars. To illustrate the trade-off between increas- 
ing R\n while simultaneously decreasing Nj, in Figure [4] we 
show the fractional uncertainty as a function of relative size 
for both Boo II and ComBer when the background density 
and A^dwarf are derived from SDSS (upper lines). In addition, 
in Figure [5] we show the fractional uncertainty in luminosity, 
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which roughly correlates with 3 x a Rxp We stress that while 
the measured uncertainties are shown as function of the satel- 
lite's relative size, the total number of stars is also chang- 
ing. In particular, Nj varies from higher values on the left to 
lower values on the right. This figure shows that better results 
are obtained for larger FOV, with the decrease in uncertainty 
slowly plateauing as the FOV increases. This improvement is 
a consequence of the increase in the total number of stars (see 
also §5). In p ractice howe v er, for regions much larger than 
those used by Mar tin et al.l (|2008) the improvements can be 
further limited, due to potential background fluctuations and 
gradients across the field. We therefore conclude t hat it is not 
possib le to improve substantially the results of iMartin et al.1 
(2008) by choosing a smaller FOV. The ultimate limiting fac- 
tor for a SDSS-like survey is the very low density contrast 
of the satellites with respect to the background, a reflection 
of both their low surface brightness and the relatively shal- 
low limit of SDSS, a limitation that can be applied also to 
SkyMapper, expected to reach a similar depth than SDSS. 

In summary, for SDSS-like surveys, our simulations sug- 
gest that large FOVs will yield more accurate results, as long 
as the background density remains homogeneous across the 
field. 

3.2. LSST 

The Large Synoptic Survey Telescope (LSST) will image 
the full sky every two nights to a single exposure depth of 
r = 24 over a period of ten years, with a fi nal co-added mag- 
nitude limit of r ~ 27 (llvezic et al.ll2008l) . Because LSST is 
expected to find UFD satellites out to much larger distances 
as compared to SDSS, the same issues will apply in regard 
to measuring structural parameters in the limit of low num- 
bers of member stars. An added challenge at LSST magnitude 
limits will be unresolved background galaxies which will in- 
creasingly contribute to the bac kground source counts fainter 
than r > 22 (e.g., see Fig. 1 of Reid & Maj ewskilll993l) . To 
assess the expected performance of LSST we carry out mock 
observations of the UFD galaxy ComBer, placed at two dif- 
ferent distances: its true distance of 44kpc and farther away 
at400kpc. 

We estimate the number of UFD stars, AWarf, detectable 
down to the full co-added LSST survey magnitude limit using 
a theoretical luminosity function typical of the UFDs so far 
discovered: a 1 3 Gyr old population w ith overall metallicity 
of_[Fe/H]= -2.27 dGirardiet al.l 120041) . assuming a IChabrieri 
(2001) initial mass function. For ComBer at 44kpc, we ob- 
tain A^dwarf = 3420. We combine this information with the 
known angular size of this system to determine its central den- 
sity. To derive a density contrast we then need to estimate the 
background density. For a survey like LSST, this variable is 
ultimately tied to our ability to characterize the contamina- 
tion, due either to foreground Milky Way stars or unresolved 
background galaxies. We estimate the expected contamina- 
tion from Milky Way stars at these magnitudes by generating 
fake photometric data fo r our Galaxy using the TRILEGAL 
tool dGirardi et al.l 120051) . For the relevant magnitude range 
(22 < r < 27) this yields foreground star densities that are neg- 
ligible compared to those of background galaxies within our 
Color-Magnitude-Diagram (CMD) window^ (as we illustrate 
in Figure Hfb)). Our primary contamination will then come 
from unresolved background galaxies. To estimate their con- 

4 We define a region around the main sequence of the isochrone described 
above, located at 44kpc. 



tribution to T, b , we use the Hubble-Ultra-Deep-Field (HUDF) 
catalog under the assumption that the distribution of galaxies 
is isotropic. Using a region of the CMD similar to the one 
used to select UFD stars, we estimate a background density 
of Eg a i = 14.4 gal arcmin" 2 down to the limit of r = 27. This 
yields So = 2.2/(1 — //), where [i is the fraction of background 
galaxies that can be potentially removed photometrically. If 
we assume, for instance, that we can remove 80% of this 
contaminatiorQ we obtain £rj ~ 1 1 . This number may seem 
su rprisingly low com pared to the density contrast achieved 
by iMunoz et al.l (1201 Oh using shallower data. However, this 
is only a reflection of the fact that unresolved galaxy counts 
grow significantly faster at faint magnitudes. This is true even 
if we allow the density of background galaxies to vary by a 
factor a few with respect to the HUDF. 

A benefit of deeper photometry is the increase in the to- 
tal number of detections Nj. For ComBer at 44kpc, and 
assuming we can remove 80% of the galaxy counts, we get 
N T = 5050 for a half-light radius to FOV ratio of 0.25, yield- 
ing <tr 1/2 ~ 4%. The lower line in Figure |4] shows the depen- 
dence of the relative uncertainty with the satellite's relative 
size for this particular choice of density contrast value. Both 
the density contrast and the total number of stars are large 
enough that the measured uncertainty is rather insensitive to 
the choice of FOV. 

One of the unique aspects of LSST related to UFDs is the 
expectation of finding galaxies as faint as My = -1 out to 
the virial radius of o ur Galaxy, located roughly at 400 kpc 
(Tollerudet al. 2008). In Figure |6j a) we show the required 
fraction of background galaxies that must be removed as a 
function of r-band depth to achieve a given density contrast in 
the context of observing a ComBer-like satellite (My = -3.8, 
R\/i = 75 pc) located at 400 kpc. The figure shows curves for 
three different values of £o, 15, 20 and 50. To achieve a den- 
sity contrast of at least 20 for magnitudes fainter than r = 25, 
more than 90% of the unresolved galaxies must be subtracted. 
We note that while the fraction of galaxies that need to be 
removed does not increase at deeper magnitudes, the actual 
background density is a strong function of depth (Figure|6jb)). 

Alternatively, one can ask what is the limiting distance at 
which we can characterize the half-light radii of a given UFD 
with an intrinsic uncertainty of 10% or better (or equivalently, 
reach a density contrast of 20 or higher). The answer to this 
question is ultimately a function of our ability to separate ef- 
ficiently stellar from non-stellar detections. For illustration 
purposes, we will assume an upper limit of 90% of galaxy 
removafl With this number we find that for a ComBer-like 
satellite, we obtain a density contrast of 20 or more out to 
about 400 kpc, but not farther. In contrast, for fainter UFDs, 
i.e. with My ~ -2, comparable to systems like Segue I, Boo II 
or Willman I, the limit distance is 200 -250 kpc. The same is 
true for brighter (My ~ -5) but larger UFDs (i.e. systems 
with half-light radii larger than lOOpc). Achieving density 
contrasts that allow proper determination of structural param- 
eters for distant UFDs will thus present a major observational 
challenge. 

5 This assumption is based on the results of Mufioz et al. ( 2010) where we 
found that using the DAOPHOT morphological parameters \ and sharp, we 
can remove up to 80% of background galaxies down to a magnitude limit of 
r ~ 24.5. In practice, this fraction is likely to be lower for r = 27. 

6 This is in essence an arbitrary number at this time, since we do not know 
with certainty what will be possible to achieve with LSST, but we deem it 
a sufficiently challenging upper limit for the purpose of setting a limiting 
distance. 
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FIG. 6. — Upper panel: required fraction of background galaxies that must 
be removed as a function of r-band depth to achieve three different density 
contrast values: 5, 15, 20 and 50. This particular set of curves represents 
simulated data for a co-added LSST-like survey of a satellite with absolute 
magnitude My = —3.8, half-light radius of R1/2 = 75 pc located at 400 kpc. 
Lower panel: The upper line shows the number density of filtered back- 
ground galaxies used in panel (a) as a function of r magnitude. Galaxy counts 
were obtained from the HUDF. As a comparison, the lower line shows count 
densities for Milky Way stars calculated using the TRILEGAL tool. 

4. TARGETED OBSERVATIONS OF UFDS 

Follow-up targeted imaging of UFDs can provide more ac- 
curate structural parameters than the discovery survey data 
(ISand et al.ll2009t iMufioz et all 120101) . We next explore how 
to optimize targeted observations to provide the most accurate 
structural parameter estimates. In contrast to survey data, the 
practical observing constraint is the total exposure time which 
translates into a trade-off between photometric depth versus 
area covered. We choose to analyze a ComBer-like satellite, 
although the conclusions are applicable to the majority of the 
UFDs, for two scenarios requiring the same exposure time: 
a single lxl deg 2 pointing down to r = 24.5 and four point- 
ings down to r = 23 .0. This choice of area was motivated by 
Muno zet al.l (120101) who used the telescope/instrument com- 
bination of CFHT/MegaCam. 

The FOV of the instrument is now fixed, and thus R l / 2 is 
determined by the size of the satellite. For observations of 
ComBer down to r~ 24.5, the extrapolated number of stars 
belonging to the dwarf galaxy corresponds to A^warf ~ 750. 
Assuming a backgr ound density similar to that estimated by 
iMunoz et al] (12010b , we obtain N T ~ 1500, Ry 2 = 0.1 and a 
density contrast of So ~ 35. Thus, for a single deep pointing 
covering 1 deg we obtain ctr 1/2 ~ 5%. In comparison, obser- 
vations down to r = 23 covering twice the size of the FOV 
(four times the surveyed area), result in Eo ~ 10, R\/ 2 = 0.05 
and <tr ]/2 ~ 11%. This example shows that, in general, given 
the choice of depth versus area, the former yields better uncer- 
tainties. In light of this result, it is relevant to assess whether 
there is an optimal depth for which uncertainties will be un- 
der 10%. Figure [7] shows that, for a given total number of 
stars and satellite's relative size, the overall uncertainties only 
marginally improve for density contrasts beyond 20. Thus, 
in planning targeted observations of UFD galaxies, exposure 
times should be calculated so that Eq ~ 20. Deeper photom- 
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FlG. 7. — Fractional uncertainty cr Ri/2 as a function of the satellite's relative 
size, R\/2- Different panels correspond to different choices of the total num- 
ber of stars Nr. In each panel, three different density contrasts are shown. 
This figure shows that (i) density contrasts higher than 20 do not result in 
a significant improvement in the uncertainty of the measured half-light radii 
and (ii) the optimal range FOV size is such that 0.1 < R\/2 < 0.3. 

etry, presumably resulting in greater density contrast, will not 
have an appreciable effect on the overall quality of the mea- 
sured parameters. Of course, this does not preclude the ob- 
server from seeking higher density contrasts for different pur- 
poses. In addition, one should plan the area coverage so that, 

Total number of 
10-20% within 



ideally, 0.1 < R 1/2 < 0.3 and N T > 1000. 



stars as low as Nt = 500 will yield <jR ir 
the same range of satellite's relative size and density contrast 
value (panel (a) of Figure|7]i. 

5. GENERAL TRENDS 

The main goal of our simulations is to provide observers 
with a set of tools to help select the appropriate set of parame- 
ters when designing observations of UFDs. As shown in pre- 
vious sections, the distribution of recovered Ri/ 2 values will 
depend on all three parameters simultaneously: total number 
of stars, satellite's relative size and density contrast. How- 
ever, it is possible to distill these results into a set of "rules of 
thumb" by isolating the effects of each variable independently. 

In Figure Oa) we explore the effects of varying the satel- 
lite's relative size. In this figure, we show the distribution 
of R\/2 for three cases where we have varied the size of the 
satellite. For this particular example we fix Nj = 2000 while 
simultaneously keeping the density contrast fixed at So = 15. 
The general trend illustrated by panel (a) is such that, for a 
constant number of stars, in the limit where the size of the 
satellite is comparable to the FOV (^1/2 ~ 0.5) or when the 
field is too large (^1/2 < 0.1), our ability to recover the intrin- 
sic value of the half-light radii decreases significantly. This is 
readily understood since for large satellite sizes, with respect 
to the FOV, we are not able to probe the full profile, while for 
small sizes, the number of UFD stars decreases (for a fixed 
Nj) and the satellite's signal is diluted. In practice however, 
as one uses a larger FOV, the number of stars also increases, 
and therefore the measured uncertainty behaves as shown in 
FigureHl that is, no improvement or worsening in accuracy is 
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FIG. 8. — Distribution of measured R^ 2 for simulations where we vary one of the studied variables. Left panel: we vary relative sizes. The lowest value of 
R\j2 = 0.05 corresponds to a satellite occupying 10% of the FOV, while the highest value of 0.50 corresponds to a satellite covering the full FOV. Middle panel: 
we vary the total number of stars. The density contrast and satellite's size are kept constant. Right panel: we vary the density contrast. The total number of stars 
has been kept constant at Nj = 1000. and R\n = 0.2 

achieved beyond R\i 2 ~ 0.1. Observing larger FOVs should 
not be a problem, since one can always use only a subset of 
the data. 

We next investigate the effect of varying the total number of 
stars, Nt- Figure[3fb) shows the distribution of recovered Ru 2 
for three simulations with input values Nt = 500, 1000 and 
5000, bracketing values obtained for various current and fu- 
ture surveys (see Table[TJ, and constant So = 5 and Ry 2 = 0.2. 
The figure shows that when all other variables are constant, 
increasing the total number of stars improves the parameter 
estimation. For the entire grid of simulations we find that for 
a wide range of So and Ry 2 values, uncertainties better than 
10% require a minimum or Nt = 1000. 

Finally, we investigate the effects of varying the density 
contrast So, on the recovered half-light radii. In Figure [HJc), 
we show the distribution of measured half-light radii for three 
different choices of central density (in units of the background 
density), 5, 10, and 50 for satellites with Ru 2 = 0.2 and 
Nt = 1000. Figure [8] shows that, as intuitively expected, as 
the stellar density contrast between the satellite and the back- 
ground increases, the distribution of measured half-light radii 
becomes narrower and better defined (i.e., more Gaussian) 
yielding a median value closer to the intrinsic value. We find 
that for a given combination of Nt and R l / 2 , the improvement 
in the measured uncertainties is most significant for So < 20. 
For higher values of the density contrast, ctr^ decreases only 
marginally. 

6. SUMMARY AND CONCLUSIONS 

With the discovery of the UFDs we have entered a new 
luminosity regime dominated by shot noise and where tra- 
ditional methods to estimate structural parameters no longer 
yield reliable results. In this article, we have carried out mock 
observations of a suite of simulated UFD galaxies aimed at ex- 
ploring how different observational choices impact our ability 
to extract reliable structural information for these systems. In 
particular, we have investigated the effects of the total num- 
ber of stars in a given survey, central-to-background density 
contrast and relative size of the satellite with respect to the 
FOV. 

Our mock observations allow us to propose a set of simple 
"rules of thumb" designed to help plan studies of UFDs with 



existing and/or upcoming surveys as well as plan follow-up 
observations: (a) the ratio of the radius of the satellite relative 
to the FOV must be smaller than 0.3, (b) the exposure time 
should be calculated so that density contrast is ~ 20. Higher 
numbers will result in marginal gains in parameter uncertain- 
ties and (c) for observations of UFDs with photometric sur- 
veys the density contrast will be fixed by the survey parame- 
ters, and therefore the primary observational choice is the area 
around the satellite. This should be chosen so that, ideally, 
N T > 1000. 

It is worth keeping in mind that these rules of thumb are 
derived for the best case scenario, i.e., the choice of density 
profile faithfully describes the real data, the satellite is per- 
fectly circular, and the background density is homogeneous 
throughout the FOV. In practice, all three assumptions are 
likely to be broken in varying degrees, which will lead to fur- 
ther uncertainties not quantified by this study. This leaves 
room for many possible extensions to the present study. Fu- 
ture work should include, for instance, an exploration of the 
effects of recovering structural parameters using a density law 
different from that used to construct the simulations. Includ- 
ing intrinsic elongation for the simulated satellites is also war- 
ranted given the high ellipticity found for most UFDs, as it is 
extending the analysis to satellites that present small tidal dis- 
turbances. 

We have shown that the new regime of extreme luminosi- 
ties probed by the UFD population often results in unreliable 
structural parameters, especially when relatively shallow pho- 
tometric surveys, such as SDSS, are used. In the particular 
case of the half-light radii, on which mass estimates depend 
linearly, numbers based on SDSS catalogs can have intrin- 
sic uncertainties as high as 100% for the faintest and more 
diffuse UFDs. These higher uncertainties will, for instance, 
propagate directly into the predicted 7-ray fl uxes coming 
from UFDs due to dark matter annihilation (e.g. Striga riet al.l 
2008). Assessing the reliability of the measured structural pa- 
rameters must be regarded as critical in these studies. 
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